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With an increasing complexity of nanoscopic systems and the modeUng thereof, new theoretical 
tools are needed for a reliable calculation of complex systems with strong electronic correlations. To 
this end, we propose a new approach based on the recently introduced dynamical vertex approxi- 
mation. We demonstrate its reliability already on the one-particle vertex (i.e., dynamical mean field 
theory) level by comparison with the exact solution. Modeling a quantum point contact with 110 
atoms, we show that the contact becomes insulating already before entering the tunneling regime 
due to a local Mott-Hubbard transition occurring on the atoms which form the point contact. 
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Introduction. In recent years, we have seen a tremen- 
dous experimental progress in the direction of man-made 
nanostructures. For example, with the wizardry of mod- 
ern semiconductor technology quantum effects in quan- 
tum dots could be revealed [T]-[7] ; in the area of molecular 
electronics transport through single molecules can now be 
studied [8r;10j ; and for magnetic storage technology nan- 
oclusters of transition metals on surfaces become relevant 
^llj . In all three examples electronic correlations play a 
decisive role since the restriction to nanostructures brings 
the electrons close to each other so that their mutual 
Coulomb interaction becomes large (compared to their 
kinetic energy or tunneling rates). As a matter of fact, 
electronic correlations are not only genuine to nanostruc- 
tures, but they also make them fascinating, both from 
the basic research point of view, with new physics oc- 
curring, and from the point of applications since strong 
correlations result in spectacular physical properties. An 
example is the Kondoesque physics which overcomes the 
Coulomb blockade and which has been observed in the 
conductance of quantum dots [2]-[5l [7] as well as for small 
clusters and individual adatoms p!2UT5] . 

The theoretical modeling of strong correlations in 
nanostructures attached to some environment (bath) 
such as the source and drain electrode in case of the quan- 
tum dot or the surface for the transition metal cluster is 
hitherto based on generalizations of the Anderson impu- 
rity model f5| l7l [TBHIS] . However, if one is not only deal- 
ing with a single or two " sites" (say the number of quan- 
tum dots), the numerical effort to solve the correspond- 
ing Anderson impurity model becomes prohibitively ex- 
pensive. More precisely, the effort grows exponentially 
with the number of "sites" for an exact [19l |20] or nu- 
merical renormalization group [18^ treatment. This re- 
stricts these methods in effect to 0{2) sites coupled 
to a bath. Related dynamical matrix renormalization 
group (DMRG) approaches [2T1 might allow for slightly 
larger systems but ultimately suffer from the same non- 
polynomial problem, except for truly one dimensional 
geometries. Potentially more efficient quantum Monte- 



Carlo methods [22 on the other hand exhibit a growing 
sign problem with increasing system size. Hence a good 
theory for correlated nano-systems with even a few cou- 
pled nano-objects is presently missing. Such a theory 
is however mandatory since future technological applica- 
tions will require the engineering of complex networks of 
such nano-objects - be it for a quantum computer or for 
a von Neuman computer based on such small structures 
so that quantum effects are no longer negligible. 

On the other side, dynamical mean field theory 
(DMFT) [23H26] along with its cluster [22] and diagram- 
matic extensions such as the dynamical vertex approx- 
imation (DrA)[27] and the dual fermion approach [28] 
has been applied to strongly correlated electron systems 
with great success: On the model level, among others, 
the Mott-Hubbard transition [24l |25], magnetism [29] . 
and kinks in strongly correlated systems [3QH32] could be 
better understood or have even been discovered. Merg- 
ing DMFT with density functional theory in the local 
density approximation (LDA) [33ll36] turned out to be a 
breakthrough for the calculation of actual materials with 
strong correlations. By construction, these DMFT cal- 
culations are done in the thermodynamic limit, i.e., for a 
macroscopically large system. 

There has been one attempt by S. Florens [37] to es- 
tablish a nanoscopic version of DMFT. The idea behind 
Florens' nanoDMFT approach is the DMFTesque limit 
of a large number of neighbors with a central site in the 
middle, surrounded by many neighbors which in turn are 
coupled to many neighbors etc. In such a geometry (e.g., 
of a Cayley-type tree), one gets a recursive method where 
the inner sites depends on their outer neighbors but not 
vice versa. Experimentally however, such a geometry 
with more and more neighbors is hardly realizable, and 
the approach has been scarcely used in practice [38 . 

Here, we hence take another route based on the DFA 
concept of the locality of the fully irreducible n parti- 
cle vertex F. While the calculations in this first paper 
will be for n = l, i.e. on a DMFTish level, we call the ap- 
proach nanoDFA in following - also to distinguish it from 
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FIG. 1: Abstract scheme of the nano-systems investigated. 
The individual nano-objects (blue sites) are connected by tun- 
neling channels with amplitude t. These amplitude might be 
rather complex and dependent on the distance as indicated 
by the t' and t" tunneling channels exemplar ily displayed. 
The whole nano system is connected via some additional tun- 
neling channels V to the environment sketched schematically 
in red/orange and assumed to be non-interacting. Possi- 
ble experimental realizations are coupled quantum dots or 
molecules connected to source and drain electrodes or a clus- 
ter of Fe atoms with the environment being a surface instead 
of the indicated geometry. Electronic correlations are par- 
ticularly strong because of the nanostructuring. Lower panel: 
Geometry of the quantum point contact with 110 atoms stud- 
ied below. 

the aforementioned nanoDMFT. Below, we introduce the 
approach, validate its range of applicability against the 
exact numerical solution for system sizes where this is 
still possible, and demonstrate its potential by hands of 
calculations for a quantum point contact with 110 sites. 

Method. As pointed out in the introduction, we are 
interested in a nanoscopic system consisting of nano- 
objects (sites) i which are hybridized via t^j, interacting 
by a Coulomb repulsion Ui and coupled by Vi^k to some 
non-interacting environment, see Fig. [l] The Hamilto- 
nian hence reads 

ija i 

+ Vi^,kclj^^^ + h.C. + eukllkaKka^ (1) 

ivkcT vkcT 

where (c-^) and /J^^ {Kk'[) denote the creation (anni- 
hilation) operators for an electron with spin a on site 
i and in lead v state k with energy e^y/c, respectively. 
While we consider a single band situation in the follow- 



FIG. 2: NanoDFA algorithm consisting of two self-consistency 
steps: 1) From the local Green function Gu and Coulomb in- 
teraction Ui the n-particle fully irreducible vertex V^i has to 
be calculated. 2) From F-^ a new Green function Gu is deter- 
mined through the Dyson and parquet equation, respectively. 



ing, Hamiltonian ([T]) can easily be generalized to include 
orbital realism, leading to an additional orbital index in 
the second quantization operators and orbital matrices 
in the Green functions and self energies below. 

As an exact solution of Hamiltonian ([T]) is possible at 
most for 0(10) interacting sites, we here propose an ap- 
proximate DFA solution. To this end, we first need to 
calculate the fully irreducible n-particle vertex on ev- 
ery site i with the interacting Green function Gii{uj) 
and Coulomb interaction Ui as an input, see Fig. [2] In 
practice, this is done by numerically calculating the cor- 
responding n-particle vertex of an associated Anderson 
impurity model. Note that the effort for this computa- 
tionally most expensive step only grows linearly with the 
number of sites and is easily parallelizable. From the 
n-particle vertex in turn, we recalculate the Green func- 
tion and proceed with the first step until convergence. 
In the case n = 1 , the one-particle fully irreducible vertex 
is simply the self energy Y{uj) which is directly related 
to Gii(uj) through the Dyson equation given in Fig. [2] in 
matrix notation for the site indices ij. For n = 2, one 
needs instead to use the parquet equations to go from 
the irreducible vertex to the reducible one and the exact 
equation of motion to get the self energy before proceed- 
ing with the Dyson eq., similarly as discussed in Refi27l 
for an infinite system. 

Let us note that the approach becomes exact in several 
limits: (i) t/^0, (ii) F^oo, (iii) number of connections 
to neighbor sites Z^oo, and, if each site couples to its 
own lead (iv) t ^ 0. While the exact quantum Monte 
Carlo simulation is impossible for large clusters due to 
the so-called sign problem, our method is sign problem- 
free. For n = 2 the approach also takes into account the 
Cooperon diagrams so that weak localization physics is 
explicitly included, as are spin fluctuations. 

Validation vs. exact result. As a first test case and to 
validate the approach for n = 1 against the numerically 
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FIG. 3: (Color online) On-site spectral function "A(0)" vs. 
hybridization V, comparing nanoDFA (lines) with the exact 
solution (symbols) for two different geometries: hopping to 
two nearest neighbors only and hopping to all neighbors. In- 
set: Conductance between two opposite leads of the benzene 
ring I.e., between the dark red leads of Fig. [l] for both topolo- 
gies of neighborhood. Together with the low off-slte-dlagonal 
self energy, the results show that the nanoDFA approach Is re- 
liable If either the hybridization strength Is sufficiently strong 
or If the Inter-cluster hopping Is to sufficiently many neigh- 
bors. 



exact solution, we consider the 6-slte benzene geometry 
of Fig. [l] (upper panel) with a constant density of states 
p In the contacts from —DioD around the Fermi level 
{D = 2t^ t = 1 sets our unit of energy In the follow- 
ing); and a site-diagonal hybridization Vijyk = V6ijy. Two 
topologies are considered: (1) hopping tij restricted to the 
two nearest neighbors In the hexagon "(nn t)" and (11) an 
equal hopping to all sites "(all t)". Both DFA and the 
exact solution are calculated by means of Monte-Carlo 
simulations [41] for U = 5t, temperature T = O.Obt. 

Fig. |3] shows the calculated local spectral function 
In an Interval 2T around the Fermi level (set to zero) 
"A(0)" = / dcjA(a;)/[e^/2T^g-a;/2T]^ rj.^^ results clearly 

show that nanoDFA Is reliable both for a hybridization 
V and If enough neighbors (In our case 5 In the "all f 
topology) are Involved In the hopping processes. There 
are some deviations for the nearest-nelghbor-hopplng- 
only case ifV since In this situation non-local correla- 
tions such as those Involved In forming a two-site singlet 
are relevant. If the hopping Is to all neighbors though, 
no deviations could be Identified all the way down to the 
last calculable point V = O.St. Below V = O.St^ the sign 
problem becomes too severe [average sign (9(10~^)], and 
only the nanoDFA solution Is possible - In a situation 
very favorable for nanoDFA because of the many neigh- 
bors. Let us note that another Indication for the small 
correction Is the off-diagonal self-energy, which Is smaller 
than 10~^t for V > t and nearest neighbor hopping as 
well as for all V and hopping to all neighbors. 
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FIG. 4: (Color online) Conductance across the quantum point 
contact of Fig. [l]( lower panel) vs. hybridization Iqcp between 
the two atoms of the quantum point contact. The conduc- 
tance becomes almost zero for |tQCp|^0.5t. Inset: Spectral 
function layer- by- layer across the quantum point contact for 
two tgcp's below and above the conductance Increase, re- 
vealing a local Mott-Hubbard "transition" In the two atoms 
forming the QPC. 

In the Inset to Fig. |3) we present the conductance G 
through the benzene-like nanostructure, depicted In the 
upper panel of Fig. [l] from one side of the molecule to the 
opposite side, calculated along the lines of Ref.42, As the 
local spectral function, the conductance again shows the 
reliability of the calculation already on the one-partlcle 
vertex level, with discernible deviations from the exact 
solution only for a low number of neighbors (I.e., 2 "nn 
only") and small hybridization to the leads. The results 
can be understood as follows: At small hybridization 
V, we have a conductance through two tunneling lead- 
benzene contacts, leading to an Increase of the conduc- 
tance (I.e., squared tunneling rate F = 7rF^p). In 
the large V region on the other hand, a Kondo resonance 
between the Individual sites of the benzene molecule and 
the respective lead forms, which suppresses the com- 
peting Inter-benzene hopping and hence the transport 
through the molecule. If hopping between all benzene 
sites Is allowed ("all f ), this effect Is less pronounced 
since there Is a direct hopping channel between the op- 
posite sites to which the voltage has been applied. 

Quantum point contact (QPC). To demonstrate the 
suitability of the approach for more complex nanosys- 
tems, let us now consider a QPC. Experimentally It 
can be realized e.g. through a mechanically controllable 
break junction of a conducting wire, see e.g. Ref. [43l [45l 
The assumed geometry Is based on a body-centered cu- 
bic basic structure narrowed to a double-cone-llke junc- 
tions as shown In Fig [l] (lower panel). For the mo- 
ment we assume a single band which might be realiz- 
able In more complex wires such as cuprate and cobal- 
tate wires, but realistic calculations, e.g.. In the spirit 
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of LDA+DMFT j33ti36]. are certainly possible. The pa- 
rameters are: = 4t, t = -1, f -0.4, V 0.4, and 
T = 0.05t. Each calculation takes about 10 hours with 
a mildly parallelization on 25 Nehalem Intel processors 
(X5550, 2.66GHz), showing that much bigger calcula- 
tions or calculations with orbital realism are possible. 
When slowly breaking up the junction the hybridiza- 
tion (tunneling) tqpc between the two atoms forming 
the point contact will change most strongly and is hence 
varied. 

Surprisingly, we observe a dramatic reduction of the 
conductance for I^qpcI ^ 0.5t in Fig. [4j Breaking up 
the junction triggers a local Mott-Hubbard "transition" 
(more precisely a crossover) of the two atoms forming 
the quantum point, see Fig. |4] inset. Therefore the con- 
ductance drop with increasing distance between the two 
atoms is faster than the exponential decay of tgcp with 
distance. Our findings might explain similar experimen- 
tal observations [44 in transition metals point contacts 
with partially filled d-shells, where the electrons are ac- 
tually similarly strongly correlated as in our calculation. 
This effect could not have been revealed with hitherto 
employed methods [45^ such as LDA, Landauer formal- 
ism or Coulomb blockade calculations. 
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